
  #include "stdio.h"
  #include "12lplq.c"
  main()
  { int i;
    double x[6];
    double a[3][6]={ {1.0,2.0,7.0,1.0,0.0,0.0},
     {1.0,4.0,13.0,0.0,1.0,0.0},{0.0,2.0,8.0,0.0,0.0,1.0}};
    double b[3]={10.0,18.0,13.0};
    double c[6]={-4.0,-9.0,-26.0,0.0,0.0,0.0};
    i=lplq(a,b,c,3,3,x);
    printf("\n");
    if (i!=0)
      { for (i=0; i<=5; i++)
          printf("x(%2d)=%13.5e\n",i,x[i]);
        printf("\n");
      }
  }

